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ABSTRACT 

Discrepancies have emerged between the predictions of standard cold dark matter 
(CDM) theory and observations of clustering on sub-galactic scales. Warm dark mat- 
ter (WDM) is a simple modification of CDM in which the dark matter particles have 
initial velocities due either to their having decoupled as thermal relics, or having been 
formed via non-equilibrium decay. We investigate the nonlinear gravitational cluster- 
ing of WDM with a high resolution N-body code, and identify a number of distinctive 
observational signatures. Relative to CDM, halo concentrations and core densities are 
lowered, core radii are increased, and large halos emerge with far fewer low mass satel- 
lites. The number of small halos is suppressed, and those present are formed by 'top 
down' fragmentation of caustics, as part of a 'cosmic web' connecting massive halos. 
Few small halos form outside this web. If we identify small halos with dwarf galaxies, 
their number, spatial distribution, and formation epoch appear in better agreement 
with the observations for WDM than they are for CDM. 

Subject headings: cosmology:theory — cosmology:dark matter — galaxies: formation 
— galaxies: halos — methods: numerical 



1. Introduction 

Standard cold dark matter (CDM) theory appears to be in conflict with observations of clus- 
tering on sub- galactic scales. First, the theory predicts a very large number of low mass halos. 
For example, of order 500 satellites of mass greater than 10 s Mq are expected within a halo as 
massive as that of the Milky Way whereas only 11 candidate dwarf galaxies with comparable mass 
are observed (Moore et al. 1999b). In the whole Local group, out to ~ 1.5 Mpc there are only 36 
known dwarf galaxies (Mateo 1998), where several thousand would be expected in a straightforward 
interpretation of the theory. The second problem with CDM is that simulations produce halos with 
more concentrated cores than those inferred from the rotation curves of real galaxies (Moore et al. 
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1999a; Ghigna 1999). There are a number of related but more complex problems. For example, the 
predicted satellites would lead to a thickening of the disks of galaxies which is not observed. And 
the dense cores, which form very early, would cause baryonic cooling to be too efficient, leading to 
disks with an order of magnitude too little angular momentum (Sommer-Larsen & Dolgov 1999; 
Navarro & Steinmetz 2000). 

These problems may well have complex astrophysical solutions. Feedback processes such as 
heating and supernovae winds might inhibit star formation in low mass halos (Bullock, Kravtsov, 
& Weinberg 2000; Binney, Gerhard, & Silk 2000), although it is not clear why such processes 
would suppress the number of satellites without destroying them all. No convincing mechanism 
has been devised to explain the fact that most observed dwarf systems have characteristic densities 
which are considerably lower than those of dark matter halos of the same mass predicted from 
ACDM theory, and it certainly seems implausible that feedback processes would preferentially 
inhibit star formation in the denser halos. Further, dwarf spheroidals are believed to be completely 
dominated by dark matter, but their stellar profiles show no sign of a central density cusp (Irwin 
& Hatzidimitriou 1995). 

The observational situation regarding galaxy cores has been rapidly evolving. It was recently 
claimed that beam smearing effects were significant, bringing observations of the inner rotation 
curves of low surface brightness galaxies closer to the predictions of CDM (van den Bosch & 
Swaters 2000; Swaters 2000) (although there remained a discrepancy with the p oc r~ L5 central 
profile indicated by high resolution CDM simulations (Moore et al. 1999a; Swaters 2000; Klypin et 
al. 2000)). But very recently de Blok et al. (de Blok et al. 2001) reanalised the same data and added 
considerably more new data, only to reach the opposite conclusion that beam smearing effects were 
not important and the data are seriously discrepant even with the milder p oc r _1 behavior of the 
NFW profile (Navarro, Frenk, & White 1997). 

There is another potential problem with CDM, noted some time ago by Peebles (2000, 2001). 
The puzzle is that faint dwarf galaxies as measured in our vicinity (for a picture of the local 
volume see Peebles (1993) p. 26 and Schombert, Pildis, & Eder (1997)) do not fill the voids in 
the bright galaxy distribution. Rather they tend to follow an even sharper distribution than the 
bright galaxies, tracing out surprisingly thin sheets around the voids. This effect has not as far as 
we know been systematically quantified, but is interesting in connection with the satellite problem 
of CDM. That is, observationally, dwarf galaxies seem both more strongly correlated with bright 
galaxies than in CDM, whilst at the same time being suppressed within large galaxy halos. In CDM 
theory both features are hard to understand, since the smaller halos formed earlier and should not 
be strongly correlated with later forming, larger scale structure. They should fill both the voids 
and massive structures alike, as we shall verify in our ACDM simulations presented below. Even if 
supernova winds or heating processes inhibit star formation in dwarf halos, and thereby suppress 
dwarf galaxies, it is hard to see why these processes would act preferentially in the voids. Another 
clue of potential significance is that observations (Metcalfe et al. 2000) point to the formation 
of dwarf galaxies at 3 < Z < 4, later than bright galaxies. This is not what one expects in a 
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hierarchical theory like CDM where small objects should form first. 

If these problems with CDM are real, they represent a remarkable opportunity. The observed 
pattern of gravitational clustering may be revealing the physical properties of the dark matter. If 
so, this will be an invaluable clue to physics beyond the standard model. Deciphering this clue 
represents an exciting challenge in which both more refined observations and numerical simulations 
will be needed. Even if the conclusion of this work is that CDM theory survives, it will have been 
strengthened in the process. Alternative theories are valuable foils against which the successes of 
the better theory may be judged. 

The small scale problems mentioned above do not negate the remarkable successes of the 
VLm ~ 0.3, VL\ ~ 0.7 ACDM cosmology on larger scales (Bahcall et al. 1999). These include the 
abundances and observed evolutionary properties of Lyman a clouds, large galaxies, clusters and 
of course the normalization and shape of the power spectrum of cosmic microwave anisotropies. 
Therefore it seems sensible to seek a modest modification of the scenario whose sole effect is to 
damp small scale structure. 

Free streaming due to thermal motion of particles, or Landau damping, is the simplest known 
mechanism for smearing out small scale structure. This process was central to the hot dark matter 
(massive neutrino) scenario for structure formation, which failed because the particle speeds were 
too high and erased perturbations up to tens of Megaparsec scales. Structure formation was 'top 
down', with galaxies forming only through fragmentation of pancakes, and at redshifts too low to 
be compatible with observation. 

Warm dark matter (WDM) is just hot dark matter cooled down. We shall review below some 
of the many ideas in the literature as to how WDM might be produced. A particularly interesting 
possibility is that if the reheating temperature following inflation is low, then standard fi or r 
neutrinos can be warm dark matter. Historically, WDM was considered (and quickly rejected) 
as a means of reconciling computed and observed cluster abundances while maintaining a critical 
density for the dark matter (Colombi, Dodelson, & Widrow 1996). That problem is resolved in the 
ACDM model by simply lowering the dark matter density. But the new problems faced by CDM 
motivate reconsideration of WDM (reincarnated now as AWDM) with a view to damping structure 
on sub- Megaparsec scales (Sommer-Larsen &i Dolgov 1999; Hogan <fe Dalcanton 2000), requiring 
much lower particle velocities. 

It is conventional to discuss WDM in terms of a canonical candidate, a light fermion with 
two spin degrees of freedom like the neutrino, assumed to decouple adiabatically and whilst still 
relativistic. The number of relativistic degrees of freedom at decoupling determines the effective 
temperature of the warm particles today, and the dark matter density determines the WDM particle 
mass mx- We shall discuss WDM in terms of this mass, but it should be remembered that what 
really matters is the particle streaming speed, given for the canonical particle in equation (A3) of 
the Appendix. 

We discuss linear perturbation theory relevant to warm particles in the Appendix, where we 
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derive the following formula for the smoothing scale, defined as the comoving half-wavelength of 
the mode for which the linear perturbation amplitude is suppressed by two: 

The scaling given here is different from that given in the previous literature on WDM (which 
mostly refers to Bond & Szalay (1984) and Bardeen et al. (1986)). The scaling usually given does 
not properly take into account streaming prior to matter-radiation equality. 

The outline of this paper is as follows. In Section 2 we review particle physics mechanisms 
through which WDM may be formed. In Section 3, we discuss the lower bound imposed on mj by 
the requirement of early structure formation and in particular of reionizing the universe (the Gunn- 
Peterson constraint). In Section 4 we discuss the phase space density (Tremaine-Gunn) constraint, 
and show that the initial thermal velocities of the particles are only relevant to the inner cores 
of halos, scales of hundreds of parsecs. The following sections are devoted to numerical results, 
revealing several potential observational signatures of warm dark matter. 

We find that replacing cold with warm dark matter has the following effects: 

1. Smoothing of massive halo cores, lowering core densities and increasing core radii. 

2. Lowering greatly the characteristic density of low mass halos. 

3. Reduction of the overall number of low mass halos. 

4. Suppression of the number of low mass satellite halos in high mass halos. 

5. Formation of low mass halos almost solely within caustic pancakes or ribbons connecting 
larger halos in a 'cosmic web'. Voids in this web are almost empty of small halos, in contrast 
to the situation in CDM theory. 

6. Late formation (Z < 4) of low mass halos, in a 'top down' process. 

7. Suppression of halo formation at high redshift (Z > 5), and increased evolution of halos at 
lower redshifts relative to CDM. 

The first four findings indicate that WDM holds some promise as a solution to the satellite 
and core density problems of CDM. The fifth is interesting as it may solve the problem raised 
by Peebles. The sixth and seventh items point to observational tests at high redshift, where 
there is the possibility of actually seeing the pancake formation as it first occurs. Recent HST 
observations probing galaxies at high redshift seem to indicate the formation of dwarf galaxies 
after bright galaxies (Metcalfe et al. 2000), which is an argument in favor of WDM. It is perhaps 
worth emphasising that WDM succeeds both in suppressing the number of satellites in large halos, 
whilst producing a distribution of dwarf halos in a cosmic web connecting higher mass galaxies. All 
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of these effects need to be better quantified in both observational data and simulations, so that 
precise statistical tests are possible. 

It is a challenging problem at the limit of current numerical codes to simultaneously represent 
the modes responsible for structure formation (tens of Megaparsecs at least) whilst resolving objects 
as small as dwarf galaxies < 1O 8 M0. We chose the compromise of studying large boxes of 20h~ 1 
Mpc with 256 3 particles in which we compare ACDM with AWDM, with WDM particle masses of 
mx = 175eV and 350 eV corresponding to smoothing scales Rs of 2.3 and 1.0 h~ x Mpc respectively. 
These values for mx are almost certainly excluded by the Gunn-Peterson limit which we detail 
below, or by the limit from Lyman a clouds, pointing to mx > 750 eV (Narayanan et al. 2000). 
Nevertheless they are helpful for illustrating the qualitative differences between WDM and CDM 
clustering on large scales. We also performed a simulation with a more realistic value of mx = 
1.5 keV, corresponding to a smoothing scale of 0.19 Mpc, in smaller 3ft. -1 Mpc boxes with 
128 3 particles. Finite size effects are important in this box, since long wavelength modes in the 
power spectrum are not properly represented. Nevertheless the simulation shows that the effects 
found at low mx persist at higher values. A second limitation is that with the current code, 
computation time becomes prohibitive for the 256 3 boxes at redshifts below unity. Therefore our 
largest simulations end at Z = 1. The signatures we identify deserve more detailed study, and 
larger simulations will be carried out in future work. 1 

Simulations have also been recently carried out by Avila-Reese et al. (2001), Colin, Avila- 
Reese, k Valenzuela (2001), and Colin, Avila-Reese, & Valenzuela (2000). They found that WDM 
results in a reduction of the number of satellites in Milky- Way sized galactic halos, and in satellite 
halos with less centrally concentrated density profiles. The work presented here is in agreement 
with these findings. Moore et al. (1999a) simulated a ~ 1O 12 M halo and found little difference in 
the density profile between warm and cold dark matter. White &; Croft (2000) used simulations to 
examine the nonlinear power spectrum in WDM models. 



2. Warm Dark Matter 

In WDM, free streaming smears out perturbations on scales smaller than the comoving scale 
Ax over which a typical particle has traveled. In a flat FRW universe, the momentum p of a particle 
scales inversely with the scale factor a. This may be understood as a stretching of the de Broglie 
wavelength with a. The peculiar speed is given by the usual formula v p = pj \jp 2 + m\ , with 
p = mxcinr/a, where a nr is the scale factor when the particle becomes non-relativistic (i.e. when 
v p = c/y/2). The motion in comoving coordinates is given by dx/dt = v„/a ~ ca nr a~ 2 at late times. 

2 

In the matter era, where aaP, the comoving displacement Ax converges. A crude estimate of the 



1 Movies of the simulations were found to be very instructive in showing how structure formation in WDM proceeds. 
We have made these, higher resolution versions of the figures and some other results not included here available on 
the web, at http://astro.Princeton.EDU/~bode/WDM/. 
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smoothing length is obtained by multiplying the speed of the particle at matter-radiation equality 
by the comoving horizon scale at that time. Up to a numerical factor this is the smoothing scale 
usually quoted (Bond & Szalay 1984; Bardeen et al. 1986); 



However the situation is actually more complicated since for the warm particle speeds of interest the 
particles become nonrelativistic long before matter-radiation equality. There is thus a substantial 
period of time in the radiation era during which the Jeans mass is low and linear perturbations 
grow logarithmically with time, whilst the comoving displacement is also growing logarithmically 
(since a oc in the radiation era). The linear evolution of perturbations during this epoch may 
be understood using the nonrelativistic Gilbert equation discussed in the Appendix, which leads to 
the improved scaling formula given in (1) above, and which we used in the simulations described 
below. But for a crude discussion, as we give now, the the scaling of the smoothing length (2) is 
certainly adequate. 

A review of the mechanisms through which WDM may be formed and their relation to cos- 
mology can be found in Colombi, Dodelson, & Widrow (1996). Warm dark matter is most simply 
discussed in comparison with the usual hot dark matter scenario. In the hot dark matter scenario 
(with a single massive neutrino) the particle mass would be 12 eV. These neutrinos would have 
thermal speeds of v s /c ~ 0.2 at equality, and their motion would erase structure on comoving scales 
of ~ 25 Mpc. If we wish to preserve the successes of CDM theory on comoving scales greater than 
say 0.1 Mpc, we need to reduce the streaming speed by a factor of 250. (In this discussion we shall 
choose cosmological parameters VLdm = 0.3 and h = 0.65, unless otherwise stated.) 

For a thermal relic particle which is relativistic at decoupling but nonrelativistic by equality, 
the typical streaming speed at equality is v s /c oc Tx/mx, where Tx is an effective temperature 
varying inversely with the scale factor. Neglecting the subtleties of logarithmic growth in the 
radiation era, this speed determines the streaming length. To lower v s by 250 we need to either 
lower Tx or raise mj. However, if the particle decouples with thermal abundance, then we have 
px oc T x mx which must be held fixed to explain the dark matter today. Thus to obtain the desired 
reduction in v s we must lower Tx by ~ 4 and increase mx by ~ 4 3 = 64, bringing it close to 1 keV. 

The factor of ~ 4 reduction in Tx could be produced in several ways. The simplest is that 
the X particles decoupled when there were 64 times the number of relativistic degrees of freedom 
than when neutrinos decoupled, i.e. 64 x 10.75 = 688. This requires several times more degrees 
of freedom than present in the standard model (106.75) or supersymmetric versions (~ 200). A 
substantially larger number is certainly possible in theories with larger gauge groups and repre- 
sentations, especially those with extra dimensions. As far as we know, no detailed computations 
of entropy release have been performed in the latter case, but at temperatures above 1 TeV the 
large number of Kaluza-Klein modes would surely become important. If the dark matter particle 
possesses only gravitational strength couplings and therefore decouples early, it would be left with 
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a substantially lower effective temperature than the radiation fluid. A plausible candidate would 
be the gravitino (Kawasaki, Sugiyama, & Yanagida 1997). 

In fact, all that is really required is that there be some release of entropy after X decoupling. 
One possibility is that the universe became dominated by some form of matter with P < |p after 
X decoupling. For example, a mild first order phase transition with some supercooling and latent 
heat release, a rolling scalar field which dominated briefly before being converted to radiation, or 
a massive species which dominated before decaying into radiation but not X's. This must happen 
when the temperature was greater than ~1 MeV if it is not to disturb nucleosynthesis. 

Another apparently viable scenario, discussed recently (Kawasaki, Kohri, <fe Sugiyama 2000; 
Giudice, Kolb, & Riotto 2000) is that the reheating temperature after inflation Trh was low, so 
that the maximum temperature of the radiation era was in fact only a few MeV. In this case, 
if the neutrinos are only produced via electron positron annihilation, their number density can 
be considerably lower than in the usual equilibrium scenario. As Kawasaki, Kohri, & Sugiyama 
(2000) point out, although this reduces the expansion rate of the universe which decreases the 
helium abundance, there is a compensating effect because there are fewer neutrinos to convert 
neutrons to protons. According to Giudice, Kolb, & Riotto (2000), in such a scenario, ordinary 
neutrinos may be resurrected as warm dark matter. Giudice, Kolb, & Riotto (2000) give Q u h 2 = 
(m^/70keV)(TR///MeV) 3 for the density of /i and r neutrinos where m„ is the neutrino mass. 
(Unfortunately Kawasaki, Kohri, & Sugiyama (2000) and Giudice, Kolb, &: Riotto (2000) are 
apparently inconsistent in their neutrino densities by over an order of magnitude. We understand 
the origin of the discrepancy is still unclear (Giudice 2000)). Recalling that the desired momentum 
(in units where c = 1) is ~ 10 _3 m^ at equal density, it is ~ 10 3 m 1/ TRH /MeV when the neutrinos 
are produced. The neutrinos are radiated from electron-positron and neutrino-antineutrino pairs 
in the hot plasma with thermal momenta p ~ STrh assuming Trh >> m u- These equations imply 
m u ~ 3 keV, and for Q, v ~ .3 and h ~ 0.65 we obtain Trh ~ 2 MeV. The numerical coincidence 
here is very interesting, and deserves further study. Note also that in this scenario, the number 
densities of fi and r neutrinos are suppressed, so that matter-radiation equality is now at a higher 
redshift (by about thirty per cent). This is helpful in that there is more time for structure to form, 
alleviating the reionization bound discussed below. Interestingly, it also shifts the Doppler peak of 
the CMB anisotropy to lower I. 

A different alternative is that X particles carry a conserved charge like baryon number, and an 
X asymmetry is formed as a result of non-equilibrium CP violating processes. In this case we would 
have (as for baryons) px ~ emxT x with e a dimensionless CP violating coupling. Assuming the 
dark matter particles have the same temperature as neutrinos, slowing them down by the required 
amount would require their mass be increased by 250, thus we would have mx ~ 3 keV, and the 
required dark matter density is attained for e ~ .004. 

One may also consider X particles produced via entirely non-thermal mechanisms. One idea 
(Thomas 1995) is that a supersymmetric condensate simultaneously decays into the lightest super- 
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particles, which form the dark matter, and baryons, in a ratio which explains the observed coin- 
cidence in their densities today. If the condensate mass were of order 300 times the dark matter 
particle mass, and if the decay took place at around 1 MeV, then the dark matter particles would 
have gamma factors of p/mx ~ 300 at that time, and speeds of 3 x 10 _4 c at matter-radiation 
equality, equivalent to a 1 keV warm particle. 

Other examples of non equilibrium production of warm particles are resonant production of 
sterile neutrinos via oscillation from active neutrinos (Dolgov h Hansen 2000; Abazajian, Fuller & 
Patel 2001), and particle production from the decay of massive particles or cosmic defects (Lin et 
al. 2000). 

Even from this limited discussion it is clear that there is an (over) abundance of mechanisms 
for producing warm dark matter, using modest extensions of physics beyond the standard model. 

3. Early Structure Formation in WDM 

If one uses WDM to solve CDM's small scale problems there is a window of allowed streaming 
speeds: too low and the WDM is just like CDM; too high and structure formation begins too late, 
as in HDM. A particularly stringent limit is the Gunn-Peterson bound (Gunn & Peterson 1965), 
provided by the fact that quasar spectra reveal little absorption from smoothly distributed atomic 
hydrogen. Thus the intervening regions of the universe must have been almost completely ionized. 
The highest redshift quasars (at Z = 5.5 and 5.8) have only recently been discovered (Stern, et 
al. 2000; Fan et al. 2000), and have normal spectra. There must therefore have been enough UV 
radiation released from stars or AGN to fully ionize the universe by this epoch (Gunn & Peterson 
1965). 

The fraction of the baryons which could plausibly have undergone gravitational collapse and 
star formation by a given redshift is easily computed (assuming the primordial fluctuations are 
Gaussian) from the power spectrum given in the Appendix. The fraction of the universe f n i 
attaining an overdensity greater than 5 C at some redshift is just probability that 5 > S c at any one 
point. We choose 5 C to be the linear theory amplitude corresponding to collapse to zero radius 
in the nonlinear spherical collapse model. In an VLm = 0.3 cosmology, one has 5 C « 1.675 (Eke, 
Cole, & Frenk 1996), slightly smaller than the canonical critical density value of 1.69. We compute 
the fraction /„; by integrating the power spectrum and employing the linear theory growth factor 
appropriate to this cosmology. For Z > 2 the growth factor g(Z) (<5o(l + Z))^ 1 where 5q is the 
loss of growth relative to an Om = 1 flat universe, at Z = 0. The latter is well approximated 
(Lahav 1991) by 

(3) 

The result for f n [ at redshifts 7,6,5,4,3 is plotted as the solid curves in Figure (1) as a function 
of the warm particle mass, assuming the canonical neutrino-like particle and a AWDM theory 
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normalized to COBE (<r 8 = 0.9). 

When baryons are bound into stars, nuclear burning releases on order 4 MeV per baryon. To 
ionize hydrogen, 13.6 eV is needed, thus a conservative bound is that at least 3.4 x 1CT 6 of the 
universe must have undergone gravitational collapse in order to re-ionize the universe. (Even this 
limit is not completely watertight, since black holes can release more energy per gram than stars). 
A more realistic bound is obtained by assuming that the first generation of stars was typical of 
star forming regions today: for such a population nuclear burning releases on order 4000 ionizing 
photons per baryon (Madau 2000). This leads to a bound f n \ > 2.5 x 10~ 4 . A greater nonlinear 
fraction is needed if one allows for recombination in denser clumps. Gnedin and Ostriker (Gnedin & 
Ostriker 1997) find (for ACDM) that 10 photons above 13.6eV are required for each net ionization. 
Note however that a WDM universe is very smooth at early times between the high peaks, so 
clumping is likely to be less significant than in CDM (Madau 2000). 

The fraction f n \ is very sensitive to the value of a% since the latter occurs squared in the 
exponent of the Gaussian probability distribution. In Figure 1 we show the results of increasing ag 
from 0.9 (the result for a COBE- normalized n = 1 primordial spectrum) to og = 1.0. The latter 
would be appropriate for COBE normalization with a small positive tilt, Sn ~ 0.02. At present such 
a tilt cannot be ruled out, although future CMB data may eventually allow one to do so. Figure 
1 shows that a minimal warm mass of order 400 eV (700 eV) is needed for as = 1.0 (0.9), if one 
allows an additional factor of 3 for recombination. A more detailed treatment of the reionization 
limit leads to the requirement mx > 0.75keV (Barkana, Haiman, & Ostriker 2001). 

The way in which structure forms in WDM is conveniently illustrated in Figure 2. This shows 
the rms mass fluctuation au in linear theory, computed in a spherical top hat window of varying 
radius, plotted against the comoving mass contained in the window. It is given in units of the 
linear theory growth factor g(Z), normalized to unity at Z = 0. For the cosmology of interest, and 
for Z > 2, g(Z) is well approximated by (<5 (1 + Z))' 1 , where 5 = 0.778. 

Above the WDM smoothing scale the theory behaves like CDM. But as the mass scale is 
lowered, the rms mass fluctuation tends to a constant, just because the density field is smooth on 
small scales. Gravitational collapse occurs in most of the mass in objects whose scale is of order the 
WDM smoothing scale. These objects collapse when the linear theory rms amplitude au reaches 
of order unity. From Figure 2 and the values of g(Z) given above we see that the corresponding 
redshift increases roughly from 3 to 5 to 7 as the warm particle mass is increased from 0.5 to 1.5 to 
3 keV. So in WDM one expects to see a great deal more evolution of galaxy scale objects at these 
redshifts than in CDM. 

What about the formation of dense collapsed objects like the centers of bright galaxies? In the 
simplest treatment, using the spherical collapse model, these regions have an overdensity of ~ 200 
at the moment of collapse and thereafter their overdensity scales as (1 + Z) 3 . Thus halos collapsing 
at redshifts ~ 4 could plausibly attain overdensities of a few times 10 today, and with the inclusion 
of mergers this is within reach of the observed overdensities of the inner parts (r < 10/i -1 kpc) 
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of bright galaxies as inferred from rotation curves (v ro t ~ 200 km s _1 implies an overdensity of 
~ 3 x 10 5 ). Since these regions contain a mass of ~ 10 n /i _1 Mq, and the number density of such 
galaxies is ~ 10~ 2 /i 3 Mpc -3 , for £Im ~ 0.3 they constitute around 1 per cent of the mass in the 
universe. From Figure 1, we see that for warm particle masses greater than 700 eV, and with 
erg = 0.9, more than one per cent of the dark matter collapsed by redshift 4. But as we lower the 
warm particle mass, the curves steepen and the corresponding redshift falls sharply. Thus it is clear 
that warm particle masses lower than ~ 500eV are strongly constrained by the requirement that 
dark matter halos with the observed properties are formed. 

Now consider smaller objects. Dwarf galaxies, like Draco and Ursa Minor have half masses 
of ~ 1.0 x 1O 7 M and half mass radii of ~ 500pc (van den Bergh 2000). This corresponds to 
a characteristic density of 2 x 10 16 Mq/Mpc 3 and an overdensity of 5 x 10 5 for currently favored 
cosmologies; they are believed to be dominated by dark matter. They also have very low line of 
sight velocity dispersions, a ~ 9 km s _1 (Mateo 1998). Small objects like this are only formed in 
WDM through fragmentation of caustics at late times in the vicinity of larger mass galaxies. As we 
discuss in the next section, their densities are perfectly compatible with the Gunn-Tremaine bound 
(which yields a mass limit nix > 400eV for the canonical warm candidate). But how can such cold, 
relatively dense objects ever form? The most plausible location is at the center of caustics. For cold 
planar collapse, the density of the virialized sheet diverges as the normal distance z tends to zero, 
as (Peebles 1993). If a particle turns around at coordinate z at time time t c , then its peculiar 
velocity is ~ z/t c . If the turnaround time is at Z = 4, t c ~ O.UHq 1 for AWDM with VLm = 0.3. 
Then if z is the comoving scale corresponding to the mass of a dwarf, d ~ 0.07/ (1 + Z) ~ .01 Mpc, 
and we find a collapse velocity of ~ 6.5 km s _1 , similar to what is observed. Clearly this issue 
requires further study, but the estimate is not discouraging. 

It is interesting to compare inferred dark halo overdensities over the full observational range, 
from dwarf galaxies to rich clusters of galaxies. As an example of a well studied cluster, see Rines 
et al. (2000). The mass is ~ 10 15 /i _1 Mq, and the inferred half-mass overdensity is 2.5 x 10 3 in 
this cosmology. (Here we are being conservative, in including the entire infall mass outside the 
virial radius.) The usual scaling argument applied to hierarchical structure formation (e.g. Peebles 
(1993), p. 628) predicts that the overdensity for collapsed halos phalo °c M h ^ +3 ^ 2 where n is the 
effective power law for the linear theory power spectrum {P^ oc k n ) on the relevant scale. The rms 
mass fluctuation om oc M h ^ +3 ^ 6 . Comparing ctm on the scale of clusters ~ 1.5 x 1O 15 M0 and 
dwarf galaxies ~ 2 x 10 7 M Q , for A CDM one obtains an effective index n ~ —2.2. The theory 
predicts a ratio of ~ 1400 for the characteristic dark matter densities of clusters and dwarfs of these 
scales. The observed ratio is 200, which is considerably smaller. Clearly a larger data base would 
be desirable, but this comparison points to a possible further problem with CDM. 
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4. The Phase Space Constraint 

Liouville's theorem, stating that the phase space density is conserved along particle trajecto- 
ries in a collisionless fluid, provides a fundamental constraint on the clustering of warm particles 
(Tremaine & Gunn 1979), which we review below. Hogan & Dalcanton (2000) and Sellwood (2000) 
have recently discussed the relevance of this constraint to WDM, which unlike CDM has a finite 
primordial phase space density due to the thermal particle velocities. However, by the time of halo 
collapse the thermal particle speeds are small and have little influence on halo formation except 
in the very inner-most regions of massive halos. The smoothing of the primordial power spectrum 
which occurred at very early times, and which may be treated in linear theory as we do in the 
Appendix, is the main effect. 

Most of the mass undergoes gravitational collapse at modest redshifts in WDM, Z < 5. By 
this time, the warm particle speeds are very low. From (A3) in the Appendix, the rms speed for 
a 1 keV particle at Z = 5 would be only v = 0.25 km s _1 , far lower than the velocities induced 
in gravitational collapse in any of the objects we are interested in today. Thus heating due to 

gravitational infall vastly outweighs the initial thermal energy, even allowing for the increase due 

i 

to phase space density conservation (v oc pi). 

What little remains of the primordial velocities could really only affect the far interiors of 
massive halo cores. To estimate this effect, consider the infall of a warm particle falling from initial 
radius into a pre-existing halo. The principal effect of the thermal motion is to give the particle 
angular momentum, producing a centrifugal barrier keeping the particle away from r = 0; only for 
radii inside this barrier is the structure of the halo significantly altered relative to a pure CDM 
halo. This barrier is only significant for radius r smaller than r min pa r^v^jv 

max j where v max is 

the maximal radial speed induced by falling through the halo and < v 2 ^ >= | < v 2 >= ^v 2 ms 
is the initial thermal transverse velocity dispersion. Assuming for simplicity a halo with a flat 
rotation curve and density run p = p s (r s /r) 2 , v max pa ^/8irGp s r 2 ln(erj/r m j„). We should take 
the initial radius V{ to be a little larger than the radius of the shell now virialized at r s when it 
turned around at time t, so ri ~ 2r s . If the halo has just virialized, the mean overdensity in r s is 
~ 150, which implies that p s ~ 50/3;, where the background matter density pt = l/(6-7rGt 2 ). Thus 
we find r m i n ~ fv r mst where t is the time when the particle turned around and / is of order 0.1 
divided by the square root of the logarithm. For collapse at Z = 5 (t = .09HQ 1 ) this gives a scale 
rmin = 250/ ~ 25/j" 1 pc, smaller than scales of interest. Generalizing this calculation to other halo 
profiles is straightforward. For example, assuming an NFW profile (7) and r m in << r s , one finds 
Vmax ~ \/32TrGp(r s )r 2 , so the conclusions are very similar. 

If we considered infall into a less symmetrical halo, we would still obtain / < 1. Thus it seems 
safe to conclude that for warm particle masses larger than 1 keV, thermal velocities are unimportant 
in determining the structure of halos on scales of a kiloparsec or above. They may be relevant to 
the formation of very small halos - if / ~ 1, then the mass in r m in is ~ 10 5 M Q . 

In order to distinguish between the effects of reducing small-scale power and increasing thermal 
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velocities, we have studied the collapse of massive halos numerically; an illustration is shown in 
Figure (3). We evolved small boxes (length 3 h l Mpc) of 32 3 particles, with particle masses 
of 6.86 x W 7 1i~ 1 Mq. For comparative purposes, we studied ACDM, and A WDM for mx = 350 
eV, first with both the appropriate linear power spectrum and the appropriate thermal velocity 
distribution, and then with only the appropriate power spectrum but no thermal velocities. We 
do include the particle velocities from the Zeldovich approximation in the initial conditions, with 
the simulations starting at Z = 40. Identical phases were utilized for all three simulations. In the 
Figure, we have identified a particular object for study, which has a mass of 2 x lO 11 /?- 1 ^/© (over 
3000 particles). We show a phase space (r,v) plot of all the particles within a region centered on the 
highest density point. The results are shown at redshifts 8, 1 and 0. In ACDM, the halo shows little 
evolution, having already begun to collapse at Z = 8 and relaxed to virial equilibrium well before 
redshift 1. The WDM halos both show substantial evolution between Z = 8 and 1, but not much 
between 1 and 0. They are very similar to each other at all redshifts, illustrating the unimportance 
of the initial thermal velocities on these mass scales. The value of pj < v 3 > near the center of this 
halo falls by more than four orders of magnitude over the course of the simulation. 

Let us now briefly recall the bound given by Tremaine & Gunn (1979) on the possible phase 
space density of collapsed objects. This is most relevant to the objects with highest density and 
lowest velocity dispersion, namely the dwarf spheroidal galaxies. Liouville's theorem states that 
the fine grained phase space density at the location of a fluid particle is constant in the absence of 
collisions. In dense regions, the measured phase space density can decrease due to coarse graining, 
but it cannot increase. For the canonical neutrino-like particle we are considering, the phase space 
density in the initial conditions has a maximum value of hp 3 , where hp is Planck's constant. This 
is the upper bound for the collapsed object. If one assumes that after virialization the velocity 
distribution is Maxwellian, with one dimensional dispersion a 2 , one obtains a final phase space 
density of p/mx x (2ira 2 m 2 x )^^ . Thus the limit is 



With the further assumption that the core is modeled as a King profile, with core radius r c , one 
finds for the central density p = 9a 2 / '(AnGr 2 ) (see e.g. Peebles (1993), p. 435), and thus 



Draco and Ursa Minor have core radii of ~ 200 pc and a of 9 km s _1 (van den Bergh 2000), so one 
sees that mx must be greater than about 400 eV for warm particles to dominate these galaxies. 
But this bound is weaker than others we have described above. 




(4) 



In observational units this translates to 






(6) 



-13- 



5. The Simulations 

Three simulations were carried out: ACDM and two AWDM models, one with warm particle 
mass mx=350 eV and one with 175 eV (vq = 0.12 and 0.048 km s -1 , respectively). The cosmo- 
logical parameters are = 0.3, = 0.7, and h = 0.67; all models were normalized by as = 0.9 
to be consistent with observations on large scales. Each simulation is a periodic cube of length 20 
h^ 1 Mpc containing 256 3 particles, making the particle mass 4 x 10 7 h^ 1 Mq; a Plummer softening 
length of e = 1.22/i" 1 kpc was used. The time step was set to be y^0.05e/a m , where a m is the 
magnitude of the largest acceleration experienced by one of the particles. All runs were carried out 
with the P 3 M code of Ferrell & Bertschinger (1994), in the parallel version of Frederic (1997). 

To generate the initial conditions particles were displaced from a regular grid in the standard 
Z'eldovich approximation. The phases used for the displacements and velocities were identical in 
the three simulations. The matter power spectrum for the LCDM model was computed using 
the COSMICS package (Bertschinger 1995; Ma k Bertschinger 1995; Bode & Bertschinger 1995). 
For WDM this was modified by the transfer function \T^\ 2 of equation (A8) in the Appendix, 
suppressing small scale power. Warm particles were also given additional random velocities in 
accordance with the appropriate Fermi-Dirac distribution. 

The number and distribution of virialized halos are central to our study. To locate collapsed 
objects the HOP code of Eisenstein & Hut (1998) was used, smoothing over 32 particles when 
calculating the density at each point and using an outer density contour of 100 times the mean. Each 
halo thus identified was then considered as an isolated object, and those particles not gravitationally 
bound were removed. Then, using the most bound particle as the halo center, the radius r2oo of the 
sphere enclosing a mean overdensity of 200 was found, as well as the mass M200 inside this radius. 
Objects with M200 < W 9 h~ 1 M@, i.e. with fewer than 25 bound particles, are removed from further 
consideration. 

These simulations do not have the resolution required to go to the smallest observational 
scales, since we only resolve objects > h^ 1 10 9 M Q . A second limitation is that the runs become 
computationally costly at low redshifts. Hence we shall show results for these large boxes only 
down to redshift unity. We have performed smaller (128 3 particles with 3/i -1 box size) simulations 
to redshift zero, for a more realistic warm particle mass of 1.5 keV, which we also discuss below. 

6. Evolution and Clustering Pattern 

Figure (4) shows the projected density in three runs of ACDM and AWDM with warm particle 
masses of 175 and 350 eV. From top to bottom, the boxes are shown at redshifts Z of 3, 2 and 1. 
The most obvious distinction is that the warm boxes have less structure at Z = 3, especially on 
small scales. This illustrates the constraint posed by the requirement of early galaxy and quasar 
formation in a WDM cosmology, and of course the requirement that the universe be re-ionized at 
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least by Z ~ 6. 

The second feature of the figure is that the clustering pattern on large scales is identical in the 
three cases, as is to be expected given the identical phases used to generate the initial conditions. 
However on smaller scales interesting details are revealed. In the ACDM model the regions between 
the 'cosmic web' structure are filled with small halos, which do not occur at all in the warm models. 
Small halos are formed in the warm models, but only within the caustic surfaces ('ribbons') of the 
'cosmic web'. The different clustering patterns are shown in more detail in Figure (5), where the 
spatial distribution of halos of three mass ranges are plotted. In the CDM model the small halos 
M < 9 x 10 9 /i~ 1 M Q fill the voids, whereas in the warm models they closely trace a 'cosmic web' 
pattern of remarkably thin (and cold) caustic sheets. 

An illustration of the typical environment of a halo is shown in Figure (6). Here the local 
density around each halo is calculated by measuring the total mass within a 1 Mpc sphere 
surrounding the halo center. The cumulative fraction of halos as a function of the local density is 
plotted. The plot shows that half of the less massive halos (10 9 to 9 x itfh^M®) in the cold model 
are located in low density regions at less than twice the mean density. In the warm models, almost 
all such halos are in regions denser than this. The median local overdensity increases from 1.5 for 
CDM to 4 and 7 for 350 and 175 eV WDM respectively. Results for larger halos only (> 1O 12 M ) 
are also plotted, revealing little difference in their mean environments between the three models. 
This is due in part to the fact that the halo itself is included when calculating the local density; 
centering a sphere on a massive halo ensures that the enclosed mass is a few times overdense. In 
the 175 eV model, the distribution for more massive halos is quite similar to that for small halos, 
showing that smaller halos form in the same filamentary environment as larger ones. 

A blow-up of the projected density of one octant of each cube is shown in Figure (7). The 
suppression of small scale structure in the warm models is evident. The closer view of massive 
halos reveals many satellite halos in the ACDM model, but very few in the warm models. This 
is quantified statistically in Figure (8), which shows the average number of satellite halos within 
each halo. This number defined by measuring the distances separating each pair of halos: if this 
distance is smaller than the sum of their radii (defined as the radius enclosing mean overdensity 
200), then the less massive halo is counted as a satellite of the more massive halo. There is strong 
suppression of the number of satellites of the massive halos in WDM, by roughly a factor of 5. 
Somewhat surprising is that the effect is even larger for the 350eV warm particles than for the 175 
cV particles. The number depends on both the rate at which satellites are accreted and the rate 
at which they are tidally destroyed (Colin, Avila- Reese, & Valenzuela 2001). In WDM small halos 
are found in denser regions; on the other hand, having formed later these are less dense and thus 
more easily disrupted when encountering a massive halo. 

Careful examination of the web structures in the central column reveals at late times a fairly 
regular pattern of instability with wavelength of maximum growth being comparable to the trans- 
verse ribbon dimension. We believe this fragmentation pattern is a real reflection of the top down 
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instability; this is discussed in more detail in the next section. 

To summarize the findings of this section, we see that the bulk of the low mass halos are formed 
in the classic 'top down' structure formation scenario, via the fragmentation of large ribbons and 
pancakes. However we see this taking place on relatively small scales, within a scenario in which 
hierarchical 'bottom-up' structure formation is simultaneously occurring on larger scales. In the 
next section we quantify the abundances and properties of halos in the three simulations. 

7. Halo Abundances and Structure 

Figure (9) shows the halo mass function found in the three models. There are two features 
distinguishing the AWDM models from ACDM. First, there is a strong suppression of low mass 
halos, consistent with the cutoff seen in the linear power spectrum. This suppression is seen for 
masses below the scale set by equation (1), ^irpi,Rg, or 4 x 10 11 and 4 x 10 12 h~ 1 M Q for mx =350 
and 175 eV respectively. Secondly, the mass function becomes much steeper again at the lowest 
masses. However, these low mass halos are not formed early, as the corresponding cold dark matter 
halos are. Rather they are formed via pancake fragmentation as explained above, and therefore at 
lower redshifts and with correspondingly low densities. 

Jenkins et al. (2001) provide a general formula which predicts the abundances of collapsed dark 
matter halos for a range of cosmologies (the approach of Sheth, Mo, & Tormen (1999), improving 
on the standard Press-Schecter excursion set approach, gives quite similar results). We compare 
this prediction with the numerical results in Figure (9). ACDM is fit very well across the entire 
mass range of the simulation. For the WDM models, which have a strong maximum in the power 
spectrum within the wavenumber range involved in nonlinear clustering, the excursion set approach 
fits poorly on scales below that set by the smoothing length R$. The rms mass fluctuation a on 
scales relevant to small masses is roughly constant, as is shown in Figure (2). Thus the effective 
spectral slope is highly negative, and well outside the parameter space in which the analytic formulas 
have been tested (see Figure (9) of Jenkins et al. (2001)). These formulae also do not predict an 
upturn at smaller scales, which is hardly surprising, since the formation of the small halos occurs 
via caustic fragmentation, and has very little to do with the initial power spectrum smoothed on 
the appropriate comoving mass scale. 

This last point concerning caustic fragmentation is worth examining in more detail, as it may 
be questioned how realistically this has been simulated — there may be numerical relaxation within 
caustics which leads to spurious collapse. If such relaxation is important, it could be reduced 
by increasing the softening length. As a test, the mx=175 eV model was again evolved with a 
Particle- Mesh only code; this new simulation is thus essentially collisionless. A mesh of 512 along 
a dimension was used, making the force resolution roughly the same as the initial interparticle 
spacing, and roughly a factor of 30 lower than the the resolution of our primary runs. The resulting 
mass function is displayed in Figure (9) as a dashed line; it is quite similar to the full resolution 
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case. Although there are somewhat fewer low mass halos (as would be expected), the abrupt change 
in slope of the mass function at 1O 11 /i _1 M0 is reproduced. As objects near this mass contain over 
2000 particles, it should not be too surprising that the reduction in resolution has little effect on 
such scales. The similarity seen in the mass functions (and in visual examinations of the projected 
density) of the two different resolution runs makes us confident that the fragmentation seen in 
filaments is not a numerical artifact. 

A second line of approach to ensure that the small halos are the result of caustic fragmentation 
is to examine the evolution of the number density over time. Figure (10) shows the cumulative 
mass function at redshifts 4 and l for the ACDM and mx=350 eV models. For standard ACDM 
there is relatively little change in the low mass end between these two redshifts; most of the small 
halos are already in place, so the change is less than a factor of two. On the other hand, in the 
WDM model there is an increase of roughly an order of magnitude in the number density of lower 
mass objects. The late formation epoch of these small halos also supports the concept of their 
resulting from a top down instability. 

The differences in formation history are reflected in the internal properties of halos. The 
density profile of each halo was fit to a NFW profile (Navarro, Frenk, h White 1997) 

*) = 4^4. m 

r(r + r s ) z 

which is a convenient parameterization. The core size r s is roughly the radius where the rotation 
curve of such a halo should flatten, and p(r s ) is the value of the density at that radius, which we 
term the 'characteristic' density. This was done by computing the density in radial bins of width 3 
/i _1 kpc and finding the best-fit parameters c = r2oo/r s and r s by minimizing x 2 ; if the reduced x 2 
gave less than a 90% probability for the fit then that halo was not used in the following analysis. 

Figure (11) shows the resulting concentration parameter c, r s , and p s as a function of halo 
mass for each of the models studied. For the lowest masses the differences are quite substantial. 
The WDM halos show lower densities and central concentrations than the corresponding cold dark 
matter ones. Thus the WDM halos are in better agreement with observations of nearby dwarf 
galaxies than standard CDM (van den Bosch et al. 1999). While smaller, these differences persist 
up to higher mass halos formed in the normal (bottom up) fashion. In the normal CDM model 
the central parts of high-mass halos are from accreted low-mass, high-density halos (Subramanian, 
Cen, & Ostriker 2000); in WDM such accreted objects are rare. Colin, Avila- Reese, & Valenzuela 
(2001), using a warm particle mass of 604 or 1017 eV, found that lower-mass satellites found around 
Milky- Way sized halos have concentration factors lower by a factor of 2 compared to ACDM. This 
appears to be in good agreement with Figure (11). 
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8. A Higher Warm Particle Mass mx 

The WDM simulations discussed above are for particle masses too low to reconcile with the 
limits discussed in Section I, which indicate mx > 1 keV. While they are useful for identifying 
qualitative trends, one may ask if the differences with standard CDM discussed above persist to 
higher particle masses. To test this, an additional pair of simulations were made; one was ACDM 
with the same cosmological parameters as before, and the other was A WDM with mx = 1.5 keV. 
The box size was chosen to be 3 fa _1 Mpc, with N = 128 3 particles. These runs were carried out and 
analyzed in the same manner as described in Section 5; because of the smaller number of particles 
they could be evolved to redshift zero. For such a small box, it is inevitable that there are sizable 
finite size effects, and these are visible in the orientation of the largest pancake. 

The trends seen for lower warm particle masses are also present in this 1.5 keV run. Figure (12) 
shows the positions of halos identified in the two simulations; shown are halos with bound masses 
larger than 2.7x 10 7 /i _1 M Q (25 particles). The number of low mass halos is suppressed in the WDM 
model by a factor of over three, and it can be seen that small halos trace the caustic structures 
connecting the largest ones. This is in contrast with small CDM halos, which occur in voids as well. 
The WDM halos tend to be found in denser environments as well. Using a sphere of 300/i _1 kpc to 
define the local overdensity, only 20% of the WDM halos are in regions with local overdensity less 
than unity, while 40% of the CDM halos are. The larger WDM halos have fewer satellites: only 
17 satellites are found around the ten largest WDM halos (with masses > lO lo /i _1 M0), while the 
ten largest CDM halos have 97. The WDM halo profiles are also quite distinct from CDM. Fitting 
with the NFW profile, the mean concentration parameter is lower by a factor of 2.5, and the mean 
core radius is larger by a factor of 3. The overdensity at r s has mean ~ 2 x 10 5 in ACDM with a 
few halos above 10 6 ; in WDM the mean is ~ 2 x 10 4 with ~ 10 5 as the maximum. 

There is definitely room for more detailed study, but this pair of runs indicates that the trends 
described in our larger runs do persist as the warm particle mass is increased. 



9. Conclusions 

We have identified several key observational features of WDM clustering which serve to dis- 
tinguish it from cold dark matter. The principal difference is due to the 'top down' formation of 
structure below a characteristic mass scale, which we can roughly quantify as 

Below this mass scale, objects are formed primarily by the fragmentation of pancakes and ribbons. 
They are rarer and considerably less dense than halos of the same mass in CDM. And their spatial 
distribution is very different - they are concentrated in sheets and ribbons running between the 
massive halos, an effect which has been noted for some time for dwarf galaxies in the local universe 
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(Peebles 2000, 2001). Finally, the dwarf galaxies form later than bright galaxies, contrary to the 
hierarchical picture of clustering in CDM theory. Again WDM seems to be in better agreement 
with the observations on this point (Metcalfe et al. 2000). 

This work points to the need for the following observational tests to verify or refute the 
predicted signatures, namely: 

• Observations revealing the ionization history and Lyman a cloud distribution at redshifts 
higher than 6. 

• Observations revealing the history of galaxy formation and checking whether low mass galax- 
ies formed later than higher mass galaxies, as in WDM, or earlier as in CDM. 

• Determination of the abundance, masses, density profiles and clustering pattern of dwarf 
galaxies relative to large galaxies. A recent study (Zabludoff & Mulchaey 2000) indicates that the 
dwarf to giant ratio increases with increasing density, consistent with our findings here. Likewise 
the apparent absence of dwarf systems in the voids noted by Peebles (2000, 2001) needs to be be 
carefully quantified and compared with CDM and WDM simulations. 

• High accuracy observations measuring the inner parts of galaxy rotation curves and compar- 
ing them with detailed N body simulations of AWDM and ACDM. 

• Observations of the abundances of satellites in galaxies other than the Milky way. Especially 
important is a determination of the characteristic densities of dwarf dark matter dominated systems, 
as compared to dwarf halos in simulations. 

It is clear that a great deal also needs to be done in simulating WDM, to better quantify the 
qualitative features we have noted here. This will require higher resolution codes and considerable 
amounts of computation time. 
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A. Appendix: Linear Theory 

In this appendix we review the standard results concerning thermal relic warm dark matter and 
then discuss a simplified approach to understanding damping by free streaming, which illustrates 
the physics and tells us how the power spectrum should scale as the mass of the warm particle, 
and the cosmological parameters are varied. We explain the calculation of the smoothing length 
Rs given in equation (1) above. 

For a thermal relic X which decouples when relativistic, the abundance relative to photons is 

where g* dec is the effective number of relativistic species present at X decoupling, contributing the 
entropy of the universe. Bosonic degrees of freedom contribute unity to g*, fermionic degrees of 
freedom contribute |. The first factor in (Al) represents the dilution of X particles relative to 
photons as the gdec degrees of freedom present after X decoupling decay into photons, electrons, 
positrons and neutrinos (which together have g* = -%■). The second factor is the further dilution of 
nx /n^y when electrons and positrons annihilate. Finally, the last factor is the number of degrees of 
freedom gx contributing to the number density. Bosons contribute g = 1 but fermions g = |. This 
is divided by the number of degrees of freedom for photons, which is 2. A light neutrino species 
contributes 1.5 to g, and it is conventional to use this as the fiducial value for gx for the warm 
particle. 

The abundance of X's fixes their density through px = mxnx- This yields 

g dec 1.5 keV 

As mentioned in Section 2, if entropy producing processes occur subsequent to X decoupling, they 
act in a manner as to mimic an increase in g& ec . 

The second quantity of interest is the velocity dispersion of the X particles, since this is what 
causes the smearing of small scale primordial perturbations. If the X particles decouple when 
relativistic, the distribution function (e( p//Tx ) + l) -1 remains constant until gravitational clustering 
begins. All particle momenta scale as a -1 which we can describe by scaling Tx accordingly. When 
the particles become non-relativistic we can use p = mxv, with v their speed. The distribution 
function is then (e^ ) + 1) _1 , with 

at redshift Z. The rms velocity is 3.571t>o. 

A crude estimate of the free streaming lengthscale is obtained by calculating the comoving 
distance a particle moves after matter-radiation equality. This is proportional to (A3) evaluated 
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at Z eq = 2.39 x 10 4 J7m^ 2 , times the comoving horizon scale at that time, 16.0/(f2M^ 2 ) Mpc (for 
three light neutrinos). Thus one obtains the scaling given in (2), from Bond & Szalay (1984) and 
Bardeen et al. (1986). As mentioned in the introduction however this scaling is incorrect due to 
additional streaming in the radiation era, as we shall now discuss. 

For the masses mx we discuss here, most of the smoothing of the initial perturbations occurs 
when the particles are nonrelativistic. The free streaming may then be understood via an analytic 
solution of the Gilbert equation obtained by treating the background radiation as uniform, and 
making a certain approximation to the Fermi-Dirac distribution (Brandenberger, Kaiser, & Turok 
1987). This approximation is useful for revealing the main parameter dependence of the transfer 
function. 

The perturbation in each plane wave mode of X is governed by Gilbert's equation (for a 
derivation see Brandenberger, Kaiser, & Turok (1987)) 



where ^^(z, zq) is a solution to the free streaming equation in the absence of self-gravity. F(z, z') = 
ln(l + £ /_1 ) — ln(l + z _1 ), where z = r/(4r*), and r is the conformal time. The scale factor is chosen 
to be unity at matter-radiation equality, and is given by a = 4(z + z 2 ) (Brandenberger, Kaiser, 
Sz Turok 1987) (Note that z is not the redshift). In (A4), the quantity A = kT*vo(Z eq ) where 
k is the comoving wavenumber, and the comoving value of the scale (1 + Z eq )r* = 19.4fi^ h~ 2 . 
The comoving horizon scale given earlier is 2(y/2 — 1) times this scale. Finally, vo(Z eq ) is given in 
equation (A3). 

Equation (A4) is derived by approximating the distribution function [er + l)~ x as e~r times 
an appropriate normalization factor (Brandenberger, Kaiser, & Turok 1987). In this approximation 
bosons and fermions behave identically. 

Equation (A4) shows the scale T*vo(Z eq ) entering, in comoving units. This is the speed of 
the warm particles times the horizon at matter-radiation equality. This leads to the estimate of 
R >eq in (2). However, the Gilbert equation requires initial conditions, and specification of the 
non-gravitating solution S^^z, zq). This is where additional parameter dependence enters. The 
point is that no growth of perturbations occurs whilst the X particles are relativistic, because the 
Jeans length is of order the horizon scale. When the X particles slow down, the Jeans length starts 
increasing, reaches its maximum, and thereafter decreases. The Gilbert equation only applies after 
the maximum Jeans length is attained, which is when the X's become nonrelativistic. The X 
particles' rms speed falls to ~ 0.3c when z becomes greater than 



(with rax hi keV), which we adopt as the initial value in the Gilbert equation. Note that this 
combination of parameters differs from that entering R >eq in (2). 




(A4) 



z w m^lh 2 /m x )Hl.h/g x )\ 



(A5) 
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The free streaming solution ^^(z, zq) should be calculated using the full relativistic equations, 
by integrating the perturbations across horizon crossing. However, no additional cosmological 
parameter dependence enters this phase since it is deep in the radiation era. If we are only interested 
in the suppression of structure relative to that in cold dark matter, it is reasonable to take a simple 
isotropic ansatz (i.e. a perturbation of the local temperature) for <^-°{^(z, z ) and then compare its 
evolution with that predicted from the same equation for the k = mode, which is unaffected by 
streaming. In the same approximations as used above, we find 



1 - \A 2 F(z,zp) 
(l + A 2 F(z,z ) 2 ) 3 ' 



With the source fully specified, equation (A4) is easily iterated numerically into the matter era 
where the amplitude of the growing mode density perturbation is compared to that for k = to 
obtain the transfer function. 

The smoothing kernel depends on the smoothing scale R >eq in (2), and zq in (A5). Di- 
mensional analysis then yields a smoothing scale R >eq f{x) with f(x) an arbitrary function and 
x = (VLh 2 /mj) 4 (1.5/ gx) • The results of the Gilbert equation code are well fitted by 

T k x = (1 + (akfy 5 (A7) 

with a = 0.05(Ox/-4)- 15 (V-65) L3 (keV/mx) L15 (1.5/c/x)- 29 and k in h Mpc -1 , i.e. f(x) = .05^ 
and 4/3 ~ —0.183. This additional dependence, which weakens the decrease of smoothing length as 
mx is increased, is due to the fact that higher mass particles become non-relativistic sooner and 
therefore go through a longer period of free streaming in the radiation era. 

A full Boltzmann code calculation (Ma 2000), yields a fit 

T k x = (l + (akr)- 5/U (A8) 

with 

a = 0.048 (Ox/-4)- 15 (V-65) L3 (keV/m x ) 1 - 15 (1.5/5x)- 29 , (A9) 

and v = 1.2, which is accurate to a few per cent over the relevant range of k. The scaling, and 
the fitting function, are illustrated in Figure (13). The warm dark matter power spectrum to be 
input into N body simulations is then given by that for cold dark matter multiplied by \T k x \ 2 . The 
particle velocity distribution is given in (A3). 
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Fig. 1. — The collapsed fraction f n i, denned to be the fraction of matter in the universe which 
has attained a linear theory perturbation amplitude corresponding to collapse to zero radius in 
the spherical collapse model. The fraction f n \ is plotted at redshifts 7, 6, 5, 4, 3 (lower to upper 
curves) for warm particle masses from 0.2 to 2 keV. The solid curves are for the theory normalized 
to us = 0.9 and the dashed curves for a$ = 1.0. The horizontal line shows an estimate of the 
fraction of mass which must be collapsed in order to re-ionize the universe. This is based on the 
ionizing radiation released by typical stellar populations today, multiplied by an estimated factor of 
3 to account for inefficient propagation and recombination. The observation of a quasar at Z = 5.8 
whose spectrum exhibits little absorption due to neutral hydrogen thus sets a lower limit on the 
mass of the warm particle, mx > 400 eV or 700 keV for as = 1.0 or 0.9 respectively. Cosmological 
parameters VLm = 0.3, S7a = 0.7 and h = 0.65 are assumed. 
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Fig. 2. — The rms mass fluctuation gm in linear theory, in a spherical top hat window, given in units 
of the linear theory growth factor g(Z) which is normalized to unity today. The solid lines show gm 
against the comoving mass in the window for WDM particle masses of 0.2, 0.5, 1.0, 1.5, 2.0, 3.0 and 
5.0 keV (bottom to top). The theories are normalized to Gs = 0.9, consistent with an assumed n = 1 
primordial spectrum and the COBE results. The dotted line shows the same quantity gm / g{Z) for 
CDM. Cosmological parameters = 0.3, f^A = 0.7 and h = 0.65 are assumed. 
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Fig. 3. — The evolution of the particle distributions in phase space. A small halo of mass 2 x 
lO n /i -1 M has been selected for comparative study in (left to right) ACDM, AWDM, and AWDM 
power spectrum but without thermal velocities. From bottom to top: Z = 8, 1, and 0. 
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right: ACDM, m x =350 eV and m x =l75 eV A WDM. Top to bottom: redshift Z =3, 2, and 1. A 
simulation with mx ~ 1 keV would have an appearance intermediate between the left and central 
columns. (A higher resolution version of this Figure is available at the web site referred to in the 
introduction.) 
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Fig. 5. — Position of gravitationally bound halos at redshift Z = 1, corresponding to the top row 
of Figure (4). Bottom to top: high (M 20 o > 9 x lO n /i _1 M ), intermediate (9 x lO 9 /i _1 M < 
M 20 o < 9 x lO 11 /i" 1 M ), and low mass (1O 9 /i _1 M < M 200 < 9 x 1O 9 /i _1 M ) halos. Left to right: 
ACDM , m x = 350 eV and m x = 175 A WDM. The radius of each circle is r 20 o- In WDM, the 
formation of low mass objects is in ribbon-like structures connecting higher mass halos. Note the 
almost complete absence of small halos outside these regions, while the highest mass halos (bottom 
row) are quite similar to LCDM. 
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p/p 

Fig. 6. — The mass density in the vicinity of halos. The plot shows the cumulative fraction of 
halos as a function of the mean density in a 1 h^ 1 Mpc radius sphere centered on the halo. Thick 
lines are for halos with mass (9 x 10 9 > M > 1O 9 /i _1 M ), thin ones for only mass greater than 
1O 12 /i _1 M . One sees that the WDM models have far fewer halos in low density neighborhoods: 
they form instead in higher density caustics. For high mass halos the differences in environment 
are relatively small. 
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Fig. 7. — Projected density in one octant of the 20h Mpc cubes shown above. The small scale 
structure present in the CDM model is clearly visible, as is the relative suppression of satellite 
galaxies about the larger halos in the warm models. A higher resolution version of this Figure is 
available at the web site referred to in the introduction. 
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Fig. 8. — Average number of satellites, as a function of the mass of the parent halo. See text for 
method of satellite identification. 
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Fig. 9. — Number density of halos greater than a given mass for the three models (thick lines). 
Top to bottom: ACDM, m x = 350eV (v = 0.048 km s" 1 ) and m x = 175eV (v = 0.12 km 
s _1 ) AWDM. The thin dashed lines are predictions from the formula of Jenkins et al. (2001). The 
dashed line is a repeat of the mx = 175eV run using a PM only code. 
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Fig. 10. — Cumulative mass function at redshifts Z = 4 and Z = 1. Solid lines: ACDM; dashed 
lines: m x = 350eV A WDM. 
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Fig. 11. — Mean values for the NFW concentration parameter c = r2oo/r s , core radius and core 
density (denned in equation 7) as a function of mass M200 Filled squares: ACDM. Crosses: 
m x = 350eV (v = 0.048 km s" 1 ) A WDM. Circles: m x = 175eV (v = 0.12 km s" 1 ) A WDM. 
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Fig. 12. — Position of gravitationally bound halos at redshift Z = in a 3 h~ Mpc box, for ACDM 
(left) and AWDM with mx = 1-5 keV (right). The radius of each circle is r2oo- The smallest halo 
mass is 2.7 x 10 7 /i~ 1 Mq. Finite box size effects are visible in the orientation of the largest pancake, 
but it is clear that the environment of small halos is distinct between the two scenarios. 
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Fig. 13. — The transfer function Tj? = \J / P^ DM as computed with an accurate Boltzmann 
code (Ma 2000), for m x = 0.2 0.5 0.75, 1.0 and 2.0 keV, (bold dashed curves, left to right). The 
1 keV curve is then re-plotted after rescaling k by m 1 ^ 15 , for each mass, as solid lines. The same 

4 

curve is also re-plotted after the rescaling of m|-, as dotted lines, demonstrating that this scaling 
is incorrect. The fit we have adopted in equation (A8) is shown as a long-dashed curve. 



